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Abstract 


This  paper  presents  a  perspective  on  the  initial  transient 
problem  in  steady-state  simulation.  In  particular,  it  enumerates 
five  generally  accepted  facts:  1)  Conditions  prevailing  at  the 
beginning  of  a  simulation  influence  sample  paths.  2)  The  extent 
of  influence  is  a  function  of  the  strength  of  autocorrelation. 

3)  Some  initial  conditions  are  less  detrimental  than  others  are. 

4)  Truncation  reduces  bias  but  usually  increases  variance.  5)  So 
far  no  complete  solution  exists.  The  remainder  of  the  paper  describes 
a  proposal  for  solving  the  problem.  It  relies  on  the  relatively  weak 
assumption  that  the  conditional  means  in  a  stochastic  process  of  inter¬ 
est  are  related  linearly.  An  estimator  of  the  steady-state  mean  is 
described  which  has  considerably  less  bias  than  one  can  achieve  via 
conventional  truncation.  An  interval  estimator  is  also  described  which 


follows  from  standard  regression  theory.  A  test  for  residual  bias  is 
presented  which  enables  a  user  to  judge  whether  or  not  sample  data  meet 
the  minimal  requirements  for  the  proposed  technique  to  apply.  A  second 
test  allows  a  user  to  judge  whether  or  not  a  more  efficient  estimation 


technique  can  be  used. 
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1.  Perspective 


As  you  the  audience  know,  the  chairman's  charge  to  the  Invited 
speakers  today  is  to  provide  an  up-to-date  perspective  on  the  problem 
of  the  Initial  transient  in  steady-state  simulations.  Although  each 
speaker  has  his  own  distinct  view  of  this  problem,  one  hopes  that  all 
can  agree  on  a  skeletal  characterization  of  it.  In  particular,  in  a 
simulation: 


1.  The  initial  conditions  that  prevail  at  the  beginning  of  a 
run  influence  the  sample  path  that  each  stochastic  process, 
represented  In  the  simulation,  follows. 

2.  The  extent  to  which  the  initial  conditions  affect  a 
stochastic  process  at  a  given  point  in  a  run  is  a  function 
of  the  degree  of  autocorrelation  in  the  process. 

3.  Some  intial  conditions  influence  a  stochastic  process  at  a 
given  point  in  a  run  to  a  lesser  extent  than  other  initial 
conditions  do. 

4.  Truncation  of  observations  near  the  beginning  of  a  run  re¬ 
duces  the  bias  in  the  sample  mean  as  an  estimator  of  the 
steady-state  mean  but  generally  increases  its  variance. 

5.  No  completely  satisfactory  procedure  for  resolving  the 
problem  has  appeared  yet. 

In  support  of  Point  5  note  that  during  the  past  year  three  of  the 
panelists  have  proposed  solutions.  See  Adlakha  and  Fishman  (1980), 

Kelton  (1980)  and  Schruben  (1979,  1980).  However,  I  doubt  anyone 
would  claim  a  definitive  solution. 

Perhaps,  one  of  the  reasons  that  a  solution  has  alluded  us  is 
that  while  the  published  characterization  of  the  underlying  structures  in 


»  .  ► 


-2- 

the  stochastic  processes  of  Interest  has  enabled  us  to  conceptualize 
the  problem  it  has  been  too  narrow  to  allow  a  satisfactory  solution 
for  a  wide  range  of  cases.  In  particular,  I  refer  to  the  first-order 
autoregressive  representation  used  in  Fishman  (1972)  and  more  recently 
in  Kelton  (1980)  and  Schruben  (1980).  To  overcome  this  Inadequacy  we 
describe  a  generalization  of  that  model  and  show  how  it  leads  directly 
to  an  estimator  of  the  steady-state  mean  that  is  relatively  free  of 
contamination  from  an  initial  transient.  The  remainder  of  this  section 
provides  a  concise  description  of  the  essential  points  of  the  proposal. 
The  Appendix  contains  a  more  detailed  exposition. 

2.  Proposal 

Let  I  be  the  initial  conditions  that  prevail  in  a  simulation  and 
let  X^,...,X^  be  the  sample  record  collected  on  the  ith  of  2m' 

independent  replications.  The  objective  Is  to  estimate  the  steady-state 
mean 

u  *  11m  E  X^  1=1 . m  =  2m'. 

j-~  J 

For  convenience  of  exposition  we  take 

Ni  <  Ng  1  m'. 


n21-l  *  N1  +  k 


n21  *  N2  +  * 


The  more  general  case  of  arbitrary  n<|...(n  Is  described  in  the 
Appendix. 
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We  now  Impose  a  restriction  on  {X^;  j=l ,2 .... >  that,  while 
weak,  has  relatively  profound  Implications  for  the  estimation  of  y. 

If  the  regression  of  »Xj-2, * • *  on  Xj  is  11near  and  of  the 
form 

E  (Xj  I  X1,...,Xj_1,  I)  =  b  +  a$  Xj_$  (*) 


then 


^(i)  1  _n-s  /  -s  ^ 

k  n  =  ^  Xi  p  <  k  <  n,  i=l . m 

K,ni  n1  K  j=k+1  J  1 


has  expectation 


c  A 

E  y 


(1) 

k,ni 


9hi 
n^  -k 


(**) 


where 


st  ■  Oh1'11) 


0  <  Y  <  1 ,  p  <  t. 


Here  k  is  the  truncation  parameter.  More  importantly,  the  estimator 


2  rm  ,w  a(21)  m  A(2i-1). 

y  "  ^  (N2  uk,N2+k  '  N1  ^♦k  ) 


j 


E  "  *  “  ♦  (9N,+k  -  »N2+k>  '  (N2  -  Nl> 

■  V  *  o  (YNl+k'p/(N2-N1)). 


has  expectation 
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Observe  that  whereas  the  within-replication  truncated  sample  means 
(**)  have  bias  0  [y*  ^/N^)  the  new  acroes-replication  estimator 
(***)  dilutes  this  bias  to  0  ,  provided  that 

N2  s  21^.  Moreover, 

» 

vSr  y  =  2a2  (N1  +  N2)  /  m  (Ng  -  N^2, 

a2 

where  a  is  defined  in  (19)  in  the  Appendix,  gives  an  asymptotically 
(N,-**)  unbiased  estimator  of  var  y.  Provided  that  {y!^  ;  i  =  l,...,m} 

~  nr~r  ’i 

are  normal,  one  can  treat  (y-y)  /  ^fvar  y  as  t  distributed  with  m-2 
degrees  of  freedom. 

Let  us  now  concentrate  on  the  plausibility  of  (*)  as  an  underlying 
characterization.  Clearly,  (*)  hold  for  autoregressive  processes  with 
normal  disturbances.  It  also  holds  for  a  variety  of  stationary  sequences 
of  nonnegative  random  variables  with  gamma  marginal  distributions.  See 
Lewis  (1979).  It  also  holds  for  Markov  chains.  More  generally,  provided 
that  the  expectation  on  the  left  is  bounded  for  all  j,  one  can  treat 
the  right  side  of  (*)  as  an  approximation  to  the  left  side  that  either 
becomes  exact  for  some  p  or  whose  error  of  approximation  can  be  re¬ 
stricted  by  making  p  suitably  large.  If  the  analogy  with  fitting  a 
pth-order  polynomial  is  kept  in  mind,  p  needs  to  accomodate  the  smooth¬ 
ness  requirements  of  E  (X^  I  Xj , . . . ,Xj_j ,  I)  for  j  =  p  +  1 ,p  +  2,...  . 

To  assure  oneself  of  the  relative  insignificance  of  the  bias  in 
y,  a  test  based  on  X^2jj|;  1  *  l....,m'}  is  provided  in  the 


> 
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Appendix.  Briefly,  if  E  -  E  xj|[j  s  °,  then  g^+k  *  0. 

Although  it  is  difficult  to  choose  a  k  such  that  gk  =  0,  one  would 
hope  to  be  capable  of  choosing  +  k  so  that  the  last 

observation  collected  in  each  other  the  shorter  replications,  is  unbiased. 
If  no  significance  is  found  then  y  is  the  best  linear  unbiased  estimator 
of  y  based  on  ^.l^+k*  vj^+k’  1  a  l**...m'}. 

Occasionally  one  may  pick  k  sufficiently  large  so  that 
g^/N^y  is  relatively  incidental.  In  this  case  the  estimator 


A 

U 


Y  {/]B —  *(2i-l  )  .  rn —  a(21  )  x 
i=l  ^k.N-j+k  /W2  ^k.^+k' 

m  (/N^  +  /HJ”) 


gives  smaller  variance  than  y.  The  Appendix  provides  a  method  of 
computing  a  confidence  interval  for  g^/K^y. 

Figure  1  shows  the  essential  steps  to  follow  and  procedure  M  in 
the  Appendix  contains  all  required  computational  expressions.  We  re¬ 
mark  that  the  choice  of  two  sample  sizes  N-j  +  k  and  N2  +  k  leads 
to  considerable  convenience  and  simplication  with  regard  to  estimating 
M. 


Estimate 


Collect 


'k+l’“-’  n. 
i  *  1 .  ,2m 


Check 

Significance 

°f  gN,*k 


Check 

Significance 
of  V^y 


Compute  t 

~  ~  w  a2 

Ui  gk»  and  o 


S"  Signifi- 
\  cant?. 


Use  u  and 

2o2(N1+N2)/m(N2-N1)Z 

for  interval 
estimate  ++ 


Compute 


Use  y  and 

2o2/m(^  +  sfqy 

for  interval 
_ _ estimate  tt 


Figure  1.  Estimation  of  y 


J-  See  (14)  in  Appendix  for  g^. 

++  Use  t  distribution  with  m-2  degrees  of  freedom. 
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APPENDIX 

1 .  Preliminary  Definitions  and  Assumptions 

Let  I  denote  the  conditions  that  prevail  at  the  beginning  of  a 
simulation.  Let  y.  =  E  (X .  II)  j  =  1,2,...  denote  the  conditional 

J  J 

mean  of  observation  j  and  y  =  E  X.,  the  steady-state  mean.  Assume 

J 

that 


and 


lim  y.  =  y 
j-*=°  J 


(1) 


n 

lim  l  (y.-y)  =  c 

n-»oo  j=l  J 


I  c  I  <  °°  . 


(2) 


These  relatively  weak  assumptions  are  necessary  if  we  hope  to  deduce 
a  method  of  estimating  y  when  the  influence  of  I  is  present. 

Let  us  now  assume  that  the  regression  of  X.  ,,X.  ?,...  on  X.  is 
linear  and  has  the  form 

E(Xj|X1,...,Xj_1,I)  =  b  +  as  Xj-s  (3 

j  =  p  +  1,  p  +  2,...  . 


Models  of  this  type  arise  in  conventional  autoregressive  analysis  and  in 
representation  of  sequences  of  nonnegative  dependent  random  variables. 
See  Lewis  (1979).  After  integration  over  the  domain  of  X-|,...,Xj_i  one 
has 


_r 

y .  =  b  +  )  a  y . 

J  S=1  S  vj-< 


where,  using  (1), 

b  =  y  (1  -  l  a  ). 

s=l  s 


(4) 
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Let  us  define  u>.  =  g,-g.  If  [  a„  f  1,  then  (2)  gives 
J  J  S=1  s 

a  =  -1  .  (5) 

0 

Note  the  dependence  on  the  first  p  conditional  means  g^,...,g  so  that 
c  is  a  function  of  I. 

Consider  an  estimator  of  p  of  the  form 


l  “j  l 

j=1  J 

c  =  1 — - — 


P-j 

s«0 


s=0 


1  n 

hi  ^ 


j=k+l 


p  <  k  <  n  .  (6) 


The  quantity  k  denotes  the  number  of  observations  truncated  from  the 
sample  at  the  beginning  of  the  run.  We  can  use 


"j  ■  I„,  as  “j-S 


j  >  p 


(7) 


to  characterize  the  bias  in  (6).  In  particular,  observe  that 


e  0.  =  w  + 

Hk,n  M  n-k 


(8) 


where 


9 


t 


l 


t 

j=t+l-s 


w . 


(9) 


Note  that  the  bias  In  pk  n  is  a  function  of  _p , . . .  ,pk  and 
un+l_p» • • • .un  which  are  related  to  the  beginning  and  end,  respectively. 
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of  the  sample  record  X^,...,Xn.  We  can  characterize  this  bias  more 
precisely  by  using  (7).  The  solution  to  this  pth-order  difference 
equation  has  the  form 


•  VP  J 

“j  "  ^i=1Ci  Ri  (10) 

where 

I  l^l  <  IBi+1!  <  1  1=1 . P-1 

so  that 

g  =  0(IPplt_p)  t  -  p  (11) 

and 

-  * ) = [°<"vk'p>  -  °"vni  • 

Since  g  /g.  z  0  when  n  >>  k,  we  concern  ourselves  principally  with  the 
n  k 

more  serious  bias  gk/(n-k).  In  particular,  we  now  turn  to  the  esti¬ 
mation  of  u  free  of  the  contamination  of  g^/(n-k). 

2.  Estimation  of  jj 

Consider  m  independent  replications  of  the  simulation,  each  with 

initial  conditions  I.  For  replication  i  let  n  denote  the  number  of 

A  (i ) 

observations  collected  and  p.'  '  the  sample  mean  based  on  observations 

K  , 

k  +  1  through  n^  with  k  <  n^  <  n^  j=2,...,m  and  where  the  inequality 
holds  for  at  least  one  n,.  Then 

J 


i 

I 


t » (1) 

E  UM, 


fr’rk 


1  "  1 1  •  »  •  jHl 


°m  *  I  U  l  fj 

m  i=l  j=l  3 


-1  J 

.  -  m 


Consider  the  quantities 


»-o -m'f  (f,lm 

m  .=1  i  j=1  j  K,ni 


m  m 


\  ■  «£'  iw  V  * f i>  “m, 


Then  one  has 


where 


E  9k  '  h  '  6m 


i  in  i  m  i 

A_  =  D"1  T  g_  (m  -  1  ) 

m  m  i.=1  yn.  l  ^j=1  J 

i  in  -I  m  _i 

Bm  =  Dm  l  9n  (m  “  I  ) 

m  m  t=l  1  ni  1  j=l  J 


Note  that  ^  and  Bm  contain  no  term  In  gk  ,  thus  eliminating  the  bias 
in  u  and  gk  due  to  this  term. 
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One  additional  assumption  is  needed.  Suppose  that 


n.|-»° 


Then  to  order  1/mn^  we  have 


lim  (n.-k)  var  =  a2  < 

n  ^ 


2  n-l  rm  f-l 

V  '  °  °m 


var 


var 


m 


g.  =  a2  D"1  l  f . 
k  m  i 


~  ~  2-1 
cov  (y,gk)  =  -mo  D 


m 


(17) 


(18) 


fH  in  7 

corr  (y.gj  »  -m  /  V  l  f  •  l  fl 
k  i=l  1  i=l  1 


Now,  as  n.  -*  <=°  y  and  g  converge  to  the  generalized  least-squares 
1  k 

estimators  of  y  and  g^,  respectively,  which,  in  addition  to  being 
unbiased,  have  minimal  variance  among  all  linear  estimators.  As  an 


estimator  of  a  one  has 


a2  1  vm  -  ,A(i)  r  ~  f-lx? 

°  =  55  fi  (uM(  '  “  -  9k  > 


(19) 


a2  ,2 


which  is  asymptotically  (n^  ■>  <»)  unbiased  and  for  which  (m-2)  o  /o 
has  the  chi -squared  distribution  with  m-2  degrees  of  freedom,  pro¬ 
vided  that  the  y[^  converge  to  normality.  Then  (y-y)/(o2  D"1  l  f".^)*5 
K*ni  m  1*1  1 

has  the  t  distribution  with  m-2  degrees  of  freedom,  as  does 

(5k  -  9k)/(S2  Of1  I(>1  fj. 


A 
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To  illustrate  the  estimation  method  consider  the  case  of  m  even  and 

N1  <  N2  i=l . m/2  (20) 

m  , 

l  fl1  »  m  (Ni  +  N2)/2N7N2 

Dm  »  m2  (H2  -  N1)2/4N1N2 

so  that  to  order  1 /N^ 

2 

var  \i  =  2o2  (N7  +  N2)/  m  (f<2  -  ) 

var  g^  =  2o2  N  N2  (N^  +  Ng) /  m  (N2  -  N^)2  (21) 

=  N7N2  var  |i 

7  v2 

cov  ( u ,gk )  =  -  4  N^o  /m(N2  -  N^) 

corr  (y.g^)  =  -  ^NjN2  /  (N^  +  N2)  =  /  (1+  N2/N7)  . 

Also.as  M7  -*■  «  the  distributions  of  iv-v)  lm/°  (Ni  +  N2^ 

and  (g  -  gk)(N2  -  H^)ljm/2a2  +  N?)  converge  to  the  t 

distribution  with  m-2  degrees  of  freedom.  Observe  that  a  large 
difference  N2  -  N-j •  is  most  desirable.  Later  sections  show  that  the 
sampling  plan  (20)  offers  many  conveniences  for  estimation  and  hypothesis 
testing. 


n2i-l  =  N1  +  k 


n„.  =  N  +  k 

21  2 


l  f,  -  m  (N,  +  Nj/2, 
i=l  1  1  c 


3.  Relative  Importance  of  9jc/(ni 


The  aforementioned  results  for  y  and  g^  follow  directly  from 

the  assumptions  of  a  linear  regression  in  (3),  a  conventional  form  of 

convergence  in  (17)  and,  to  a  lesser  extent,  the  asymptotic  normality 

of  yP ^  for  i  =  l,...,m.  Reflecting  on  the  method  for  a  moment,  one 
k  » n.j 

quickly  notes  that  k  does  not  serve  the  role  of  the  conventional 
truncation  parameter.  We  merely  require  k  to  exceed  p,  the  order 
of  the  difference  equation  in  (3).  The  value  of  p  depends  on  the 
smoothness  of  y.  as  a  function  of  j.  Generally,  a  choice  of,  say, 

J 

k  =  10  accomodates  a  commonly  encountered  range  of  cases  and  a  choice 

of  k  =  100  probably  includes  all  but  the  most  exceptional  cases. 

Suppose  that  for  the  selected  k  g  /( n ,  -  k)y«0.  Then  using 

k  1 


A 

y 


i=l 


(22) 


instead  of  y,  gives 


m 


var  y  =  mo  /  H  f?  ) 
i=l  1 


s  var  y  . 


To  take  advantage  of  this  opportunity  for  greater  accuracy  one  can 
test  p  =  gk/(n^-k)y  for  significance.  For  example,  suppose  we  want  to 
test  the  hypothesis  that  I  pi  does  not  exceed  a  preassigned  tolerance 
6  >0  significantly  and  that  an  investigator  is  willing  to  tolerate, 
say,  6  =  0.01  or  0.05  by  way  of  relative  error.  By  using  the 


method  of  Fieller  (1954)  for  ratio  estimates  one  can  make  the 
confidence  statement 


pr  (f^1  r,  s  p  s  f j1  r2)  =  1  ■  a 


where 


gky  -  hm  +  ( - 1 ) 1  y 

o  m  i 

;2  -  h  I  f:1 
J=1  J 


i  s  1.2 


n  in  «  a  m 

/~ 2  u  r  r-l\  /~2  u  v 


y*  =  (gk  u  ■  h  m)“  -  (p  -  h  l  fi  )  (St  "  h  I 

k  j=l  3  K  j=l  J 


h  ■  32  tz  (“>  < 

and  tm_2  (a)  is  the  1  -  a/2  quantile  of  the  t  distribution  wit! 

2 

m  -  2  degrees  of  freedom,  provided  that  y  s  0.  If  the  interval 
[fj1  r. j,  f^1  r2]  is  relatively  small,  then  one  may  want  to  use  (22) 
instead  of  (6)  to  achieve  a  smaller  variance.  For  example,  if 
|r^ |  s  6  f1  and  |r2l  s  6  f^  one  may  elect  to  use  (22). 


k  Testing  wN  A  l,  -  u  =  0 


A  remaining  issue  in  need  of  attention  is  the  implicit  assertion 
that  gn  /gk  ~  0.  To  check  the  credibility  of  the  assumption  we  test 
the  more  inclusive  hypothesis:  p  -u  =  0,  since  we  also  want  at  least 


some  observations  within  a  replication  to  be  relatively  unbiased.  At 


least  one  check  on  the  adequacy  of  this  hypothesis  seems  possible.  We 

(21-1)  (21) 

describe  the  check  for  the  case  in  (20).  Let  XN  +  k  and  +  k 
denote  the  last  observations  collected  on  replications  2i  -  1  and  2i, 


respectively  and  let  m'  =  m/2  and 


(2i  -  1)  (21) 

■  l.  —  Xu  ■ 


vi  ■  \  ♦  k 


N2  +  k 


i  =  1 ,. . . ,m'  ,  (25) 


If  N1  is  sufficiently  large  so  that  one  can  regard  +  k  -  y  =  0, 


then  Ym,  /  Vs  (Y)  /  m' ,  where 


,  m' 

v  =  _L  y  v 
m'  m'  M 


(26a) 


o  i  nr  c 

sz  (Y)  =  -j~  l  (Y  -  Y  )  , 

m'-l  i=l  1  ^ 

has  the  t  distribution  with  m'  -  1  degrees  of  freedom.  One  can 
test  the  hypothesis  at  a  prespecified  level  a  or  alternatively 
compute  the  P-value 


P1  =  Pr  <V_i*  \\'\tV  sZ  (Y)  /  m'  ) 

where  Tmi_-|  has  the  t  distribution  with  m'rl  degrees  of  freedom. 

Small  values  of  P^  lower  the  credibility  of  the  hypothesis. 

. ( 21 — 1 )  .(21)  2 

5.  Testing  N,  ver  uMi  +  k  ■  N2  w  ^  +  fc  *  o 

If  g_  /(n,  -  k)  y  is  found  to  be  significant  then  one  may  be 

111  I 

I  A/.*\  ~  2 

suspicious  of  the  assumption  that  var  w>  '  ~  o  /(n.  -  k)  i«l,...,m 

K*n^  i 


(26b) 


4 


To  test  this  hypothesis  when  the  sampling  plan  (20)  holds,  one  can 
use  the  statistic 


m'  (21-1 )  2 

F  N1  ^  K.N,  +  k  *  Ql> 

m'-l  ,m'-l  =  - jp - ^ - j- 

N2  (^k,N2  +  k  *  Q2} 


where 


n  =J_  ym'  C(21) 

Hi  ml  L  W|/  U 


S-=1  WM,  +  k 


(29a) 


Q2  =  ~W 


ym  0(21) 

Zi=1  Wk,N2  +  k 


(29b) 


If  M1  var  wj^1]  k  =  N2  var  ^k!ll2  +  k  =  ^ 

Fm»-1  m'-i  has  the  F  distribution  with  m'-l  and  m'-l  degrees 


of  freedom. 


In  the  absence  of  additional  information  one  might  incline  to 
use  a  two-tail  test.  However,  the  observation  that  the  initial  con¬ 
ditions  I  are  more  restrictive  of  the  variation  in  ^  than  they 

a(21 )  *  1 

are  of  the  variation  in  y>  *  .  L  leads  one  to  consider  a  one-tall  test. 

K *  K 

Let  F( 6)  be  the  B-quantlle  of  the  F  distribution  with  m'-l  and 
m'-l  degrees  of  freedom.  If  F(&)  s  ^  m».-|  m».-j  we  accept  the  hypothesis 
at  the  1  -  B  level. 


If  the  hypothesis  of  equality  Is  rejected  then  the  variance  ex¬ 
pression  (19)  does  not  apply.  Also,  the  use  of  the  t  distribution 


-18- 


with  m-2  degree  of  freedom  no  longer  applies  exactly.  When  (20)  is 
used  one  has 


a  N? 

u  =  f4tQ 


N2-N1  y2  N2-N1  yl 


Q,  • 


(30) 


Let 


2  a(21-1) 

°1  s  var  uk,N1  ♦  k 

2  x(2i) 

°2  =  var  yk,N2  ♦  k 


i  =  l . rtf 


so  that 


var  y  = 


m(N2-N1) 


2  ?  (N2  o2  +  N2  o|)  . 


2  2 

As  unbiased  estimates  of  Oj  and  o2  one  has,  respectively. 


.2  1  rm'  Q  »2 

*1  '  in7^  i.=]  ^k.N-j  +  k  “  V 


(31a) 


2  1  rm' 


2  m'-l 


l  (p  -  Q,)‘ 

i=l  k,N2  +  k 


(31b) 


Then  it  is  common  practice  to  treat 


Gm'-1  =  V  i  (N1  S1  +  N2  s2J 


(32) 


as  a  t  distributed  random  variable  with  m'-l  degrees  of  freedom, 
although  this  treatment  is  at  best  an  approximation. 
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Procedure  M 


Definitions:  tr(B)  =1-6/2  quantile  of  t  distribution  with  r  degrees  of 
freedom.  s  *>  tolerable  relative  bias. 


Given:  k,  N.j  <  N2,  m  (even),  m*  =  m/2,  t(n,_1(2a),  t  2(a)  and  6. 


Collect’  X{21'1}  X(2i'^  X<2i>  X<2i> 

conect.  Tyk  *Vl,...,xN2+k 


Test:  MN,+k  -  «  '  °- 


1  r»i  y"  rx*21-11  -  x(21)i 

1  •  4  IP  I  .  +Jf  AN2+kJ 


m'  1'i_1  *•  N^+k 


i  —  1 , . . .  ,m  a 


o  .2  1  vm  rv(2i-l)  v(2i)  ^2 

2*  s  m*(mr-l)  4.=1  ^XN-|+k  '  XN2+k  '  ^  * 

3.  f+?/y£s~. 

4.  If  t#  (2a),  return  indicating  failure  at  1 -a  level 

5  0(2i-D  yNl+k  x(2i-l) 

Wk’Nl+k  "l  J 


i  =  l,...,m', 


6. 


A(2i) 


N,+k 
r  2 


X 


k,N2+k  “j,k+1  "j 


(21) 


i  =  l,...,m', 


7. 


a  *  y"  s(2i-’' 

*  UM,+k  • 


a. 

B  + 

Vm'  A(2i) 

£i=1  Mk,N2+k* 

9. 

u  +■ 

(N2B  -  N1A)/m'(N2-N1) 

10. 

9k* 

N1N2  (A-Bj/m'^-Nj). 

-21- 


11.  o2  «- 


A  Cc"i(5‘-v‘  '  w)2  *  "2(°k!v«c- 


Test: 


W 


12.  h  <-  H}H2o2  [2  tm_2(a)/m(N2-N1)]2. 


2  _  2  ~2  ^(N^+Ng)  hm(N-|+N2) 


Wl 

j“lj  *  { 


is.  <  «■  Mk- 1»)‘  -  m  -  -ntif-  *  9k  - 


14.  IT  y  <  0  go  to  18  . 


ij  gk-hm-y 

15  r,  «-  -* - - - 

1  u  -hm(N1+N2)/2N1N2 


16.  r, 


v  gk-hm+y 

y2-hm(N1+N2)/2N1N2  * 


17.  If  max  (I  I,  |r2 1)  <  6  goto  23 
Relative  bias  is  significant  at  1- a  level, 

18.  v$r  (£)  -  o2  (lyN^/mCN^N,)2  . 

19.  var  (gk)  <-  2N1N2  var  (y). 

20.  Ij  ♦  V  -  tm_2  (a) Vvar(y). 

21.  I2  ♦  -Ii  +  2y. 
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22.  Return  with 

point  estimate  y, 

1  -  a  interval  estimate  [I-j.Ig]. 

bias  estimate  g^, 

variances  vSr  (y)  and  var  (gk). 

Relative  bias  is  not  significant  at  1-a  level. 


23. 

A 

V  < 

/^l  A  +  / N2  B 

H-  - —  ~ 

m' (/N  ^  +  /JT2) 

24. 

A 

var 

(y)  •<-  4o2/m  (/TJj-  +  /frp 

25. 

h 

-  W  “  V2  (ct)  Mr  (fi). 

26. 

*2 

■*-  -Ij  +  2y. 

27. 

Return  with 

point  estimate  y, 

1  -  a  Interval  estimate  [I-j.Ij,], 
variance  v£r  (y), 


indication  that  g^/N^y  =  0  at  1-a  level. 
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